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SUMMARY 

A function dependent on the solution of a set of differ- 
ential equations containing adjustable parameters, can 
be minimized by systematic search procedures in par- 
ameter space. Such procedures can be implemented 
by a hybrid system consisting of a general purpose 
analog computer and a digital expansion providing 
parallel logic building blocks. Programming of such 
a system is illustrated for a simple search procedure 
in n parameters. 

INTRODUCTION 

Thf solutions of many problems benefit from the use 
of a combination of both analog and digital computing 
devices. The required hybrid techniques of computa- 
tion can be carried out by the linkage of an analog com- 
puter with a general purpose digital computer. For 
iterative solution of problems, in particular system 
optimization, a more desirable approach is the addi- 
tion to the analog computer of an EAI HYDAC Series 
350 Digital Operations System designed to provide 
primarily: 

( a.) Flexible automatic logic control of an analog 
computation or a long sequence of analog computations. 

(b.) Large capacity, fast-access storage for prob- 
lem solutions requiring a sequence of computations 
with variable function storage. 

The assembly of many analog and digital computing 
components in an integral system gives the program- 
mer complete flexibility in his organization of a prob- 
lem solution. To demonstrate what is possible with 
a system this paper considers a problem solution which 
makes good use of the flexible automatic logic pro- 
vided by the proposed digital expansion system. For 
simplicity in explanation, this more basic feature of 
the hybrid computer is applied to the solution of a 
straightforward problem in optimization. 

The Series 350 HYDAC D. O. S. provides "logic build- 
ing blocks" which can be connected to form in this case 
the automatic control of the optimization program. 
These digital modules are electronic (solid-state) de- 
vices providing high computing speed. Any logical 
function could be implemented by relays and stepping 
switches, as has been done in many analog computer 
applications in the past. 6 ' 7 > 8 However, higher 
speeds, flexibility and programming convenience are 
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required in order to take full advantage of the speed 
of analog computers. This is not done in the repeti- 
tive operation mode where successive computations 
are simply repetitious or differ only by small per- 
turbations. The full possibilities are realized by 
speeding up the usual real-time operation where the 
operator makes intelligent decisions between compu- 
tations, on the basis of the fresh information just ob- 
tained. Successive computations can then be dis- 
tinctly different. 

The digital modules are interconnected on a remov- 
able patch panel. Programming is conveniently sim- 
ilar in principle to the approach familiar to analog 
computer users. 



Series 350 HYDAC Digital Operations System 
GENERAL FEATURES 

The Digital Operations System building blocks use two 
voltage levels to represent 0 and 1. A connection 
carrying a 1 level is said to be energized. Logical 
functions can be mechanized by interconnection of 
gates operating on these levels. Sequential logic re- 
quires flip-flops. Changes in the state of a flip-flop 
are clocked by a central system clock. Frequency of 
the clock pulses is adjustable from several megacycles 
down to manual control, pulse by pulse, by a push- 
button. This enables the operator to check parts of 
the set-up at leisure. 

Counters provided with the system can be used to gen- 
erate signals at binary or decimal submultiples of the 
central clock frequency. 

LOGIC BUILDING BLOCKS 



For the present purpose, only four types of building 
blocks are required: (See Figure 1) 

(a.) AND gates (2 inputs), the negated output is 
provided. 

(b.) OR gates (up to 6 inputs), the negated output 
is provided. (AND and OR gates, although logically 
interchangeable through negation, are distinguished 
for ease of intuitive programming). 

(c.) General purpose flip-flops (GPFF). 
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These are clocked RST flip-flops. The S (set) input 
commands the 1 state; the R (reset) input the 0 state; 
the T (trigger) input commands state reversal. Only 
one of these inputs may be energized at a time. The 
input levels are sampled by the clock pulse. Changes 
in state take place and the new state is established be- 
fore the following pulse. Assembled in groups of four, 
these flip-flops have two group- input controls: a 
"clear" input which resets all four, and an "enable" 
input which can be used to inhibit the RST inputs. 

(d.) Quadruple flip-flops in shift register connec- 
tion. (QFF) These are sets of four flip-flops with in- 
ternal connection into a shift register. Each set has 
eight outputs, the 0 and 1 states of each flip-flop. 

The inputs are as follows: 

Serial Inputs : The serial set (SS) and serial reset (SR) 
inputs are applied only to the first flip- flop and can be 
patched, in particular, from either the output of the 
last flip-flop of another QFF package for construction 
of long registers or the output of the last flip-flop in 
the register for closure of a loop. 

Shift Input (SH): Transfers the state of each flip-flop 
to the next in order on one clock pulse and enables the 
serial inputs to control the first flip-flop. 

Clear Input: (CLR) Resets all four flip-flops. 

Set and Load Inputs: One set input(S) is provided for 
each flip-flop. The load input (L) is common to the 
building block. When the load input is energized, flip- 
flops for which the set input is energized will go into 
the 1 state. 

PROGRAMMING OF THE SERIES 350 D. O. S. 

The techniques for programming the Digital Operations 
System are presently unfamiliar to most analog and 
digital computer users. The program is not intro- 
duced as a list of successive instructions but as a set 
of patching connections. This implies a task of organ- 
ization similar to that needed in the design of a digital 
computer. However, there is considerable difference 
in the emphasis placed on component economy. Econ- 
omy is essential for a computer design but of little 
concern in the programming of the digital expansion 
system as long as enough components are available. 
Thus, an empirical approach to digital expansion sys- 
tem programming, by successive simplifications and 
rearrangements is entirely acceptable. 

The logic building blocks provide entirely parallel op- 
eration though it is possible to perform digital opera- 
tions serially when desired. The digital and analog 
programs are introduced by their respective patch 
panels and interconnected by trunk lines forming one 
single set of parallel analog- digital components. 

Several signals from analog comparators can be ac- 
cepted simultaneously and combined to form several 
logical functions, while simultaneously several com- 
mands to analog switches (electronic relays, analog 
memory, mode control) can be generated. 



This paper shows some aspects of this programming 
technique, for consideration by analog computer users. 
To this effect, it considers the general optimization 
problem which is both significant and representative 
of the class of applications requiring logic operations 
rather than large capacity storage. 

For clarity, no attempt is made to reduce the number 
of building blocks to a minimum. 

USE OF SUBROUTINES 

Operations in the analog computer have to follow a 
definite time sequence. Whenever a sequence of op- 
erations has to be repeated many times, it is defined 
as a subroutine. One way to implement a subroutine 
is to use a ring counter, acting effectively as an elec- 
tronic stepping switch, where the sequence of states 
determines the sequence of operations. A ring counter 
is obtained by assembling QFF units into a shift reg- 
ister of the desired length and then connecting the out- 
puts of the last stage to the serial inputs of the first 
QFF. Initially the first flip-flop in this loop is set to 
the 1 state and the others are in the zero state. When 
the shift input is energized each clock pulse will ad- 
vance the sequence by moving the 1 state to the next 
flip-flop resetting the previous one. After all flip- 
flops have been energized successively the next shift 
returns the ringto the initial state. Each flip-flop can 
command in parallel a set of operations, including the 
initiation of other subroutines. 

The effective use of subroutines depends on two fea- 
tures: 

(a.) Appropriate control by other routines. 

(b.) Provision for delays necessary for operation 
of analog components. 

CONTROL OF INTERLOCKING SUBROUTINES 

When a routine calls, at one stage, for the execution 
of a subroutine it is often necessary to halt the advance 
of this routine until the subroutine has gone through 
its entire cycle. This interlock can be obtained by the 
formation for a given subroutine A of 3 signals. 

sA: this signal is the call for subroutine A-, it is 
formed by an OR gate which receives signals from all 
the sources of calls, generally different positions in 
one or more other routines. 

eA: this signal is generated at the time routine A 
returns to its initial state. It is produced by an AND 
gate supplied from the last state of A and the shift in- 
put of A. The eA signal is energized for one clock 
time and is called the end signal of subroutine A. If 
a subroutine should be repeated n times, a separate 
counter (analogous to index registers) can be used to 
produce the end pulse. 

eA x : this signal is produced by an OR gate sup- 
plied by eA and by the negation sX of sA. Called the 
interlock signal, eA x is fanned out to all routines which 
call subroutine A. 
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The required interlock is obtained by supplying the 
shift input of a routine by the AND combination of its 
own call signal and the interlock signals of all the rou- 
tines for which it calls. 

DELAYS 

Some routines control electronic relays, analog mem- 
ory and mode of integrators in the analog computer. 
When exercising such control the routine must incor- 
porate delays to allow sufficient time for changes in 
computation to take place so that solutions do not in- 
clude the transient response of amplifiers and the 
charging of capacitors necessary between consecutive 
operations. 

These delays can be obtained by an AND gate at the 
shift input to the routine; this gate is fed by the shift 
command and timing signals of slow rate derived from 
the carries of the central clock counter. Different 
rates can be used in different routines. 

If delays of different sizes are required in the same 
routine, the routine itself can switch control from tim- 
ing signals of one rate to those of another. Alterna- 
tively general purpose delay units (one- shot multi- 
vibrators) can be used to inhibit shifting during the re- 
quired time span. 

The delay required while the analog computer is going 
through an operate cycle, is obtained by inhibiting 
shifts of the controlling subroutine in the operate mode; 
this takes care of the cases where the run time is vari- 
able, depending on some condition in the problem. 

SIGNALS FROM DIGITAL TO ANALOG EQUIPMENT 

Signals computed by digital components and transferred 
to the analog computer are required for the following 
purposes: 

(a.) Control of the analog computer mode. 

(b.) Control of the mode of individual integrators 
or groups of integrators. 

(c.) Control of analog storage units. 

(d.) Control electronic gates used for switching of 
operations. 

For the control of the analog computer modes (OP- 
ERATE, HOLD, INITIAL CONDITION) input termina- 
tions are provided on the digital patchboard. A pulse 
applied to any input termination will trigger the cor- 
responding mode. 

The other control functions are effected by electronic 
gates which accept digital signals. The opening or 
closing of the gate is a function of the applied digital 
level. In many cases, a GPFF will buffer these com- 
mands so that the condition of the gate is maintained 
until new signals are fed to the GPFF. One GPFF can 
control several gates in synchronism. 

SIGNALS FROM ANALOG TO DIGITAL EQUIPMENT 

Particular signals produced within the analog com- 
puter are transferred to the digital equipment. These 
signals provide information of: 



(a.) The analog computer mode, and 

(b.) The state of comparators. 

The computer mode signals are ICL, HL, and OPL 
which are energized whenever the analog computer is 
in the initial condition, hold, or operate mode, re- 
spectively. 

Comparators will supply one of two levels depending 
on the sign of the algebraic sum of their analog voltage 
inputs. These levels are normalized for direct utiliza- 
tion in the digital system. 

The "Optimization Problem" 

A set of algebraic and differential equations, me- 
chanized on the analog computer, uniquely determines 
in one computation (the run) a value E, the loss- function 
or error -function to be minimized by manipulation of 
n parameters X j ..... , X n . 

This problem occurs in many situations; 

— solution of a set of algebraic equations 

— matching of boundary conditions 
--search for eigenvalues 

— process optimization 
--model building 

— curve fitting etc 

While special techniques exist for each class of prob- 
lem, the general problem is considered here. 

Analytical methods 2, 5 f or solving optimization 
problems, such as that of steepest descent, make use 
of partial derivatives determined either from the analy- 
tical problem statement or by parameter influence 
techniques'*. The experimental approach considered 
is quite different conceptually, for the search for an 
optimum is entirely based on repeated tests with dif- 
ferent parameter combinations. Experimental methods 
^> 6 can use a varied degree of sophistication to accel- 
erate convergence to an optimum and to overcome pos- 
sible difficulties due to the features of the unknown 
hypersurface. E (X^ , X n ). 

The price of refined search procedures increases very 
rapidly with the number n of parameters. For the pur- 
pose of illustrating the methods of hybrid programming 
of any search technique the following simple approach 
is selected. 

SEARCH PROCEDURE 

Assume normalization of all parameters so that a 
unique "exploration distance", h (say 1 volt on a + 100 
volt computer), is acceptable. 

Given X ^ Q and E 0 at an arbitrary starting point, the 
first partial derivatives are approximated by 

_d_E ~ » E+ i - SE" A 
d X i 2h 
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where 

*E + j = E(X l0 , . . . ,X iQ +h, . . . , X no ) -E 0 

8E- i = E(X lo ,...,X io - h ,...,X no ) -E 0 

All n partial derivatives are thus approximated by 
means of 2n calculations of E (runs). In order to 
place most of the load in the digital part of the hy- 
brid system, the values of the partial derivatives are 
stored in quantized form. Instead of storing d E/d X. , 
a quantity pj, susceptible of only 3 values, is defined 
by comparison with a limit L selected by the computer 
operator. Extension to a higher number of values poses 
no special problem. 

Pi = +1 if 8E + - SE~< - 2hL 

Pj = 0 if -2hL < 8 E + - 8E~< 2hL 

Pj = -1 if 8E + - 8E" > 2hL 

The vector with components p i defines a direction in 
n- space in which a decrease in E is probable. It is 
close to the steepest descent direction for which the 
components are dE/d\i. The program stops if all 
Pi - 0. With one of the 3 n -l possible directions de- 
termined, the parameters are changed to X io + PjAj, 
where ^ i is a fixed quantity (say 2 volts). If an im- 
provement, indicated by a decrease in E, is obtained 
the values of all parameters X^ Q are "updated" to give 
a new starting point. Additional changes by steps, 
A p in the same direction are carried out until no more 
improvement occurs. Then a smaller step-size A2 
(say . 2 volt) is selected and when this also fails, the 
partial derivatives are redetermined and the process 
repeated. If, after a new determination of the quan- 
tities Pj, even a single small step leads to no improve- 
ment, the program stops. If the analog computer does 
not overload, the program will ultimately reach one of 
the stop conditions and this should occur in the vicinity 
of the optimum. 

An advantage of this technique is that the equipment 
requirement, as demonstrated later, has a slow linear 
increase with n. It is recognized that many possible 
features of the E function (high curvatures, local min- 
ima, narrow winding "canyons", etc. ) will defeat the 
search procedure; on the other hand, the very same 
features can also defeat more elaborate and equipment- 
consuming alternatives. 

If a Euclidean metric based on voltage is used, the step 
length Ay/j p w in vary between A and&^/n depend- 
ing on the direction selected. This is acceptable be- 
cause there is no a priori reason for a voltage metric 
to have special significance. 

PROGRAM REFINEMENTS: 

The procedure outlined above can be improved at a 
small expense in equipment by addition of the following 
features: 



(1) If, after determination of direction, at least 
one large step A lis successful, then no small 
steps are used before a new determination of 
direction. However, when a single large step 
does not decrease E, a smaller step A2 is 
used. This program arrangement is a time- 
saving feature. 

(2) For non- convex E functions it is necessary to 
redetermine direction from time to time even 
if improvement is constantly obtained. To this 
end an upper limit N( say 8)isput on the num- 
ber of steps in a direction before redetermina- 
tion of direction. 

(3) The limit L for quantization of the slopes 
should be carefully chosen, If L is too large, 
the program will stop in any "flat" region of 
the E- function and for non- convex functions 
this final value of E may be far from a mini- 
mum. If L is too small, all quantities Pi will 
have values ±1 most of the time. This has the 
disadvantage of reducing the effective set of 
possible directions from 3 n -l to 2 n (for n=6 
from 728 to 64) and therefore the efficiency of 
the search procedure. 

Many refinements are possible to avoid this difficulty. 
The following is selected for its economy of equipment. 

Each slope is compared with two limit values Lj and L 2 
(L , < L 2 ). If no slope exceeds L 2 the quantities pi are 
based on L, . If any slopes exceed L 2 the quantities pi 
are based on L 2 . This procedure is poor whenever 
many slopes are clustered in a narrow range around 
either L, or L 2 , but this is unlikely to happen fre- 
quently in a long optimization. The procedure requires 
just one more comparator in the analog computer. 

(4) Whenever desired, half the exploration runs 
can be saved by replacing the slope evalua- 
tion 8E + - 8E" by 8E* , that is used forward 

2H — F 

instead of central differences. 
ANALOG PROGRAM 

From the statement of the optimization procedure the 
analog circuit (Figure 2) can be derived immediately 
as follows. 

The problem equations are mechanized by circuits 
which accept n input signals X i# The last updated 
values of the parameters Xj Q are held on n analog 
memories. The size h, A v or Ao of the increments 
SX is selected by 2 switches and distributed with both 
signs to n pairs of switches; for each parameter one 
switch determines the sign of 8Xi } the other determines 
the addition of 8^ to X iQ . This addition is performed 
by an analog memory so that the new value X j= \ iQ + 

8X^ can be transferred to the X i 0 memories when up- 
dating of the parameters is required. 



The duration of the operate period (the run) is deter- 
mined in general by the problem itself and can be de- 
pendent on the parameters. A comparator produces 
the hold signal when the appropriate condition is sat- 
isfied. At that time the value E becomes available. 
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The previous updated value E Q is held on a memory 
amplifier. The difference E - E 0 is fed to two mem- 
ory amplifiers so that ®E + and 8e " can be successively 
computed, stored and then compared. The value E 
itself is stored to permit updating of E Q when required. 
If the problem integrators can conveniently be used as 
storage for the value E the number of analog memories 
required at the output can be reduced from four to two. 

For determination of the direction parameters Pi the 
sign of 8E + - 8E ~ is sensed by a comparator, and the 



absolute value of this difference is compared with two 
constants by two more comparators. 

After a step has been carried out the change in the E 
function is stored in one of the 8E memories, say 8E + , 
and the sign of this quantity is sensed by a comparator 
to determine whether improvement has been obtained. 
Table 1 shows the conditions for digital control of the 
circuit switches. 

The comparators deliver information about the E- 
function to the digital equipment according to Table 2. 



TABLE I 



UNITS 

h/A relay 

8 X . Off/On relay 

i- 1, 2,...n 

8\i + / - relay 
Stores for E 0 and X io 
Stores for -E and - \ 

8e + store 
8E~ store 



DIGITAL LEVEL 



0 
h 

Off 



store 
track 

track 
track 



1 

A 

On 



track 
store 

store 
store 



TABLE 2 



ANALOG CONDITION 



sgn * E + 



sgn (8e~- 8e+) 

sgn ^| 8E" -8E + I -2hLj^ 



DIGITAL 
LEVEL 

0 
1 

0 
1 



NAME MEANING 



IMP 



N 



Improvement 



Best direction 
is 8X^0 



Change in X is 
advisable 
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DIGITAL OPERATIONS SYSTEM PROGRAM 

The Digital Operations System program must insure 
that, once begun, the computation progresses auto- 
matically in its search for an optimum. This part of 
the program has inputs only from the IMP, N, L^, L 2 
comparator signals and the analog computer mode 
signals. 

PROGRAM STRUCTURE 

The first step in producing such aprogram is to specify 
a flow diagram of the required sequence of operations. 
(Figure 3) This flow diagram organizes the operations 
into the following subroutines: 

Master Routine: selects the step size h, A 1 or 
A 2 tests the two stop conditions and calls alternately 
for the subroutines Explore, for determination of a 
direction, and Proceed, for carrying out steps in that 
direction. 

Explore Routine: calls for a computation of Ewith 
a single parameter displaced by +h then by -h. It steps 
through all parameters storing the j\ as they are ob- 
tained. When all pi have been obtained it sets the 8X^ 
switches accordingly. 

Run Routine: executes the computation of E for 
the prevailing V j values. 

Proceed Routine: A step size (A^ or A 2 ) being 
selected by the Master routine and a direction (a set of 
Pi) having been set by the Explore routine, the Proceed 
routine displaces the accordingly, calls for compu- 
tation of E and tests for improvement. The routine 
repeats as long as improvement is obtained, unless a 
preset number of steps is exceeded, and calls for the 
updating routine after each improvement. 

Update Routine : transfers the present values of 
\ and E into the X i 0 , E 0 memory units. 

PROGRAM MECHANIZATION 

The subroutines are instrumented in detail as follows: 

Subroutine C (Figure 4) (RUN) will start with the 
analog computer in its IC mode. The call signal sC 
is obtained by an OR gate (not shown) which combines 
all sources of calls from the EXPLORE and PROCEED 
routines. The routine will be started after new com- 
mands have been given for the 8X ; a delay is necessary 
for acquisition of the new X j by the analog computer 
(especially if some of them are initial conditions). Next 
the OP mode is triggered and the subroutine inhibited 
by the OP level. Upon automatic hold, the routine re- 
sumes with a delay to allow for acquisition of the final 
8 E values, Then the 8E+ store and, if permitted by a 
signal ENHM (enable hold minus) generated in the EX- 
PLORE routine, the 8E~ store are switched to hold by 
setting the. flip-flops controlling the memories. After 
another delay, the subroutine returns to its initial 
state, generating the end signal eC. The interlock 
signal eC x is used to inhibit the calling routines. 



Subroutine D (Figure 5) (UPDATE) will update the 
Xj 0 , E G values to the last Xi , E values. Here delays 
are provided to enable the tracking stores to acquire 
steady- state values. Since only two states are re- 
quired, the ring counter is replaced by a GPFF with 
the trigger input acting as the shift control. Signals 
eD and eD x are generated, eD* inhibits the PROCEED 
routine. 

The delays for subroutines C and Dare obtained by the 
use of a slow clock rate (SLCL) derived from a carry 
of a clock counter. 

Subroutine A (Figure 6) is the program for de- 
termination of approximate partial derivatives (EX- 
PLORE). It is associated with shift registers 8x, 8 f 
N x , N. 

8 X and N x control the analog value of 8X : The data ob- 
tained during exploration is stored in registers 8 and 
N. 

The subroutine begins with the8 x , N x cleared. First 
8i x is loaded and subroutine C called then Ni x is loaded 
(implementing 8X = -h) and subroutine C called with 
ENHM turned on. At this stage the correct value of N, 
L 1, L 2 appear on the comparators. A test is carried 
out to determine which comparison level should be 
used. The first occurence of a slope > La will clear all 
previous 8' s, set flip-flop L, which controls determina- 
tion of future 8's on the basis of L 2 alone. If Lj is not 
exceeded, 8*s are selected according to Li. The rou- 
tine shifts the new data into the 8, N stores, sends the 
8E memories to the tracking mode and repeats, shifting 
to the next parameter. When the last parameter is 
reached the 8 X , N x are cleared and the 8, Nare trans- 
ferred into 8 X , N x registers in a parallel loading op- 
eration. The end signal and the interlock signal eA x 
are generated. 

Subroutine B (Figure 7) will carry out steps in the 
chosen direction (PROCEED). It calls upon subroutine 
C for an analog run, then senses the IMP comparator 
level. If no improvement is obtained the end pulse 
eB* 1 ' is generated. If improvement is obtained, sub- 
routine D is requested for updating and the S flip-flop 
associated with the master routine is set, to record 
success. The end signal eD of the UPDATE routine 
increments a counter (r) and subroutine B repeats un- 
less r = N where the routine end pulse eB™) is gen- 
erated. The two types of end signals are combined 
to generate the interlock signal eB x . 

Master Routine M (Figure 8) first calls upon A 
for exploration, selecting step size h. Next it tests 
the stop condition 8j = . . = 8 n = 0. It then calls upon 
B (PROCEED) with the large step- size A 1# If, and 
only if, at the end of B, S = 0 (no improvement) B is 
repeated with the smaller step- size A 2 . Next the stop 
condition S = 0 is tested. If S = 1 it is reset to 0, the 
8 X , N x are cleared and routine M repeats. 



Initialization (not illustrated): It is convenient to 
initialize the digital system by clearing all its flip- flops 
from a central control. At this time the operator can 
introduce starting values for X io . The hybrid program 
will start as soon as the system clock is turned on. A 
non-repeating initialization routine is patched to per- 
form the following operations: 
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(a.) Set all flip-flops which should be energized 
initially. 

(b.) Call on the RUN subroutine to calculate and 
store the value Eo corresponding to the initial para- 
meter values. 

(c.) Start the master routine of the problem. 
EQUIPMENT REQUIREMENTS 

If no further rearrangements are carried out, the fol- 
lowing numbers of components are required, where n 
is the number of parameters. 



On the analog side: 2n + 2 switches 

2n + 4 analog memory units 
5 comparators 

On the digital side: 11 GPFF 

6 + 4 [-"-] QFF 

50 AND/OR GATES 
1 counter 



CONCLUSIONS 

The addition of theEAI Series 350 HYDAC Digital Op- 
erations System providing logic capability to the gen- 
eral purpose analog computer will enlarge its capa- 
bilities as an automatic computing device. 

The less flexible relay and stepping- switch arrange- 
ments used in the past are eliminated and full use is 
made of the computing speed of the analog computer. 
This alone brings within reach problems, such as op- 
timization, where computing times were previously 
excessive in most cases. The further addition of point 
storage in digital form on delay lines will permit the 
use of serial techniques for solution of integral and 
partial differential equations. Problems of a statisti- 
cal nature can be handled with a maximum of automatic 
operation. Many simpler auxiliary functions for 
GPAC's can be easily carried out (automatic scale 
changes, editing of graphical plots which are not con- 
tinuous in time, Fourier analysis). 
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